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1 Introduction 



The paradigm of CP violation related to the structure of CKM matrix in the Standard 
Model (and its beyond) is fueling an impressive experimental programme for the study of 
both exclusive and inclusive B-decays. The data from BaBar and Belle open new era of 
B meson physics, and open a real possibility for detecting semileptonic decay modes of B 
mesons. The semileptonic decay modes of B mesons are much cleaner samples then the 
nonleptonic decay modes since in these modes, there does not exist any problems connected 
with the presence of a third strongly interacting particle. For this reason, the study of 
semileptonic decays is one of the efficient ways for the determination of the CKM matrix 
elements. For example Vet has been determined from semileptonic meson decays [1]. CLEO 
collaboration [2] have measured the branching ratios of 7r~£~^i/ and p~ t^v 

which leads to \Vuh\ = (3.25 ± 0.14lo;29 ± 0.55) 10^"^. It is well known that for accurate 
determination of CKM matrix elements we need more reliable determination of transition 
form factors. 

In this work we calculate the transition form factors for B rjiiy decay in light cone 
QCD sum rules (LCQSR) method. The detailed description of the method and its appli- 
cations can be found in [3, 4, 5]. It should be noted that interest to the B —>■ rjip and 
B — > rj'iu semileptonic decays is due to the fact that they also give information on the 
T] — Tj' mixing (sec for example [6, 7]). 

The work is organized as follows. In Sect. 2, we derive the sum rules for the transition 
form factors ol B ^ r]£u decay. Section 3 is devoted to the numerical analysis and contain 
a summary of the results and our conclusions. 

2 Light Cone QCD Sum Rules for the B ^ r] transition 
form factors 

The B ^ Tj weak transition form factors f+{q^) and /-(g^) are defined as: 

{V{p)\uj,{l + l,)b\B{p + q)) = 2Up^ + (/+ + f_)q^ . (1) 

In this section, we calculate these form factors using light cone QCD sum rule method. 
We adopt the usual strategy of light cone sum rules method by considering the following 
correlator function: 

Il^{p,q) = tJ d^xe^'^^ (r?(^J)|^{^i(a;)7^(l + 75)6(x)6(0)^(l + 75)^^(0)}|0) (2) 

The main reason for choosing the chiral bi{l + 75)1* current instead of the bij^u current, 
which have been used in the calculation of the B ^ n weak form factor [8], is that in this 
case, twist-3 wave functions which are the main inputs in LCQSR method and which bring 
the main uncertainty to the prediction, do not contribute (see below, see also [9]). 

First let us consider the hadronic representation of the correlator. By inserting a com- 
plete set of states with the same quantum numbers as the operator bi{l + 75)?/ between 
the currents and isolating the pole term of the lowest pseudo scalar B meson, we get the 
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following hadronic representation of the correlator: 

^ (77|^7^fc|B)(g|bz75^|0) ^ {7]\u^^b\h){h\bi(l + ^,)u\0) 

= ^l(g^(p + g)>^ + ^2(g^(p + ?)2)g^ (3) 

The sum in Eq. 3 takes into account the contributions of the higher states and con- 
tinuum. Note that the intermediate states h contain not only pseudo scalar resonances 
of masses greater the m^, but also scalar resonances with — O''" corresponding to the 

operator hu. 

For the invariant amplitudes Ilj, one can write a general dispersion relation in the B 
meson momentum squared, {p + qY, as: 

^ii^f 1 [p ^ — i (hkJ. — ^ _|_ subtractions (4) 

J s — {p + qy 

where subtractions are polynomials in {p + q)^ and the spectral densities corresponding to 
Eq. (3) are given as: 

= 2U{q^)^d{s - ml) + p'As) (5) 



Pi s , , 



mil 



P2{s) = iU + f-y^^S{s - ml) + p',{s) (6) 

ITlb 

The first terms in Eqs. (5) and (6) represent the ground state B-meson contribution and 
are easily obtained from Eq. (3) using Eq. (1) and the definition 

{B\U^,d\0) = ^ 

whereas p-' represent the spectral densities of the higher resonances and the continuum. The 
spectral densities can be approximated by invoking the quark-hadron duality ansatz: 

P^is) = pf'^^eis - so) (7) 

So for the hadronic representation of the invariant amplitudes Hj we get: 

Hi = — /+(? ) sfs _|_ f — Pi — ('^) _|_ subtractions 
mb{m% - {p + q)^) J so s - (p + q)^ 

112 = 9" / — ^"^fox + / ds — — ^ \„ + subtractions (8) 

mb[ml - (p + q)^) J so s - (p + q)^ 

For obtaining sum rules for the form factors /+ and /_ we must calculate the correlator 
function in QCD. It can be done by using the light cone OPE method; i.e. by expanding 
the T product of currents near the hght cone ~ 0. After contracting the 6-quark fields 
we get 

U^{p,q)=iJ d^xe"^^r]{p)\u{x)j^{l + j,)Sb{x){l + 'y5)u{0)\0) (9) 



where Si, is the full propagator of the b quark. Its explicit expression is given by the following 
expression 



iSb{x) 



iS^{x) - igs 
1 



-ikx 



+ 



ml — k'^ 



{2{ml-k^y 



1X1/ 



Here S^ix) is the free quark propagator of the massive b quark: 



S'bix) 



ml Ki{mbV—x^] 



47r2 



i—Aj^K2{mb\/-x^) 
Att^ x^ 



(10) 



where Ki are the Bcsscl functions. 

From Eqs. (9)-(ll) it follows that, in order to calculate the theoretical part of the 
correlator function, the matrix elements of the nonlocal operators between rj meson and 
vacuum states are needed. 

Here we would like to do following remark: in all the following calculations we neglect 
the Tj — Tj' mixing since it is very small [1], and therefore we choose as the interpolating 
current for the rj meson, the SU(3) octet axial vector current: 



(12) 



In order to simplify the notation we use qV to denote :^(w7^75?i + dri',_i^/5d — 257^^755) 
and introduce notation F^, = where determined as (0|g(0)r^g(0)|?7(g)) = —ifr^qn- 

An important observation is that the terms in Sb{x) containing odd number of gamma 
matrices do not contribute, i.e. leading twist-3 terms do not give any contribution. Up to 
twist-4, the 77 meson wave functions are defined in the following way [11] : 



/•I r 1 

(?7(p)|g7M75g|0) = -ifr^p^, due-'''P''^r,{u) + —mlx'A{u) 



{vip)\'l{xhnl5gsGap{vx)q{0)\0) 



due-'^'P^ 
1 

px 



--Lml^ [ due-'^^P^'Biu) 

XocPii\ ( X^p^ 



(13) 



px 



Pa gpn 



X j Pai^x(Q;i)e-^^'^("i+""3) 



px 

(14) 



{'n{p)\q{x)gsGap{vx)-i^q{Q)\Q) = i/^mj 



P/3 Qan 



XaPft 

px 



Pa gpi. 



XfSP^ 

px 



■ipx{ai+ua3) 



X j Vai(fi±{ai)e~ 

+if,ml^{p^Xf,-pf,x^) J Vai<f\\{ai)e-'P^^^'+^^'^ 

(15) 
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where G^i/ = ^e^i.axG'^^, Vai = daida2da3S{l — ai — a2 — a^). In Eqs. (13)-(15), the 
(Pn{u) is the leading twist-2, A{u) and part of B{u) are two particle twist-4, ip\\[ai), ^p±{ai), 
<^ll(Q;i) and (pi_{ai) are three particle twist-4 wave functions. 

Here we should note that that matrix element {ri{p)\u-^^j,G°'^{ux)aapu\Q) ~ due to the 
parity invariance of strong interactions and using the identity 

Inserting Eqs. (10) and (11) into Eq. (9) and using definitions of ?7-meson wave functions 
(Eqs. (13)- (15)) for the theoretical part of the correlation function we get: 







j d^xe'""" 









327r2 



Ki{mby/-x'^) 



x^ 



1 



+ 



Va I d«e*('^^+""^>^|^FX^o(9^|| + <^||)^o(mfe%/^) 



px 



p^x^ + 2(pa;)p^ 
px 



(16) 



The next task is to carry out the fourier transformation and then to take the Borel transform 
with respect to the variable {p + g)^ in order to suppress the contributions of the higher 
states and the continuum and also in order to eliminate the subtraction terms. Finally we 
identify the same structures both at the hadronic and quark gluon levels. Subtraction of 
the continuum contribution is done by employing the quark-hadron duality as mentioned 
before. This amounts to restricting the spectral integral to s < sq. (for more details about 
this procedure see the appendix and also [8]). For the theoretical part of the correlator (2) 
function we get: 



n 



th 



2 4 
Tfl Tfl — 

"^Fv^llin + '^ll - - 2^^)p^ + V2F^^/2((^|| + (^11 + + ^A_)p^ 

F 

+2V2F^mbJ^{iprj)p^ - ^/2^mbmlJl{A)p^ - V2Fr,mlJ^{B)pi, 



V2 



[17) 



where the functions I^{lp), J^i^f), Il^{ip), and J^{^) are defined in the appendix. 

Equating the cocfiicients of the corresponding p^ and structures in Eqs. (3-6)) and 
(17) we get the following sum rules for the form factors /+ and /+ -|- /_ respectively: 



-e 



-^ilin + ^11 - - 2^^) 



m. 



U + f- = 



y^F^rrib 



o^.^ -^-V?(A) - m;j,^(S) 18) 



16 



-e M-^ 



(19) 



Eqs. (18) and (19) are our final results for the B ^ rj transition form factors. 
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3 Numerical Analysis 



Prom Eqs. (18) and (19) it follows that the main input parameters of these sum rules are 
the T] meson wave functions. The explicit expressions of the wave functions fniu), A{u), 
B{u) and </7||(Q;i), '^i.{oii), <^||(«i), and (pi_{ai) are given in [11]. 

The next input parameter of the sum rule is leptonic decay constant F^. It should be 
noted that for r] meson case, the situation is more complicated compared to the tt and 
K meson cases, due to the mixing with 77'. As we already noted that we neglect mixing 
between rj and rf . More recent analysis show that the coupling to the octet current 77 is 
= 159 MeV [13]which we use in our numerical calculations. 

In our calculations we neglect QCD radiative corrections. Therefore for consistency we 
neglect QCD corrections to the leptonic decay = 160 MeV [8] (sec also [14]). 

Having these input parameters, we can carry out our numerical calculation. According 
to QCD sum rules philosophy, first of all, a suitable range for the auxiliary Borel parameter 

should be found such that the numerical results are stable. The lower limit of 
is determined by the requirement that the terms ~ M^^" , (n > 1) remain sub dominant. 
The upper bound of is determined by requiring that the contributions of the higher 
resonances and the continuum are less then 30% of the total result. Our numerical calcu- 
lations leads that both requirements are satisfied in the region lOGeV'^ < < ISGeV"^. 
It should be noted that the light cone QCD sum rules predictions for the form factors are 
reliable at the region 



where x is a typical hadronic scale of roughly 0.5 GeV. Then q'^ < 18 GeV"^. 

In Figs. (1) and (2), the dependencies of the form factors f+{q^) and f-{q^) is 
depicted for three different values of = , 5 , 10 GeV"^ for two different values of the 
continuum threshold Sq = 35 , 40 GeV^. From these figures, we see that f+{q^) and f-{q^) 
vary weakly as varies in the region 10 GeV"^ < < 16 GeV^ up to = 18 GeV"^. 

Having the working region of M^, our next problem is to find the dependence of the 
form factors /+(/_) (g^) on g^ at given values of M^, namely at — 12GeV'^ and = 
WGeV^, at So = SbGeV^ and sq = AOGeV^. We find that /^(O) = 0.15 ± 0.3 and 
f-(0) = —0.16 ± 0.2. As we already noted, in the region g^ > ISGeV^, applicability of 
the light cone QCD sum rules is questionable. In order to extend our results to the full 
physical region, wc look some parameterization of the form factors in such a way that in 
the region, < g^ < ISGeF^, this parameterization coincide with light cone QCD sum 
rules predictions. A good parameterization of the g^ dependence can be given in terms of 
three parameters as: 



The values of these parameters for the B ^ rj transition form factors f+{f-) are obtained 
as ap = 1.07 ± 0.08 {ap = 1.1 ± 0.1) and bp = 0.19 ± 0.16 {bp = 0.28 ± 0.19) where the 
quoted errors are only due to the variations of the continuum threshold, Sq, and the Borel 
mass, M^. 



q^ <ml- 2mbX 




(20) 
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It should be noted that the form factors of B ^ r) transition can be related from 5 — > tt 
transition form factors using SU{3)f symmetry. For example, the value of f+{q^ = 0) = 
0.15 obtained via S't/(3)j? symmetry is consistent with our prediction of = 0). 

Finally, we would like to note that the background for the B — > rflv decay would be 
much smaller than that for the B — > ttIv decay, due to the much lower multiplicity, since the 
background caused by i? — r]X is one order of magnitude smaller than that of i? — > nX. 

In conclusion, we calculate the transition form factors for B rjlv decay using the light 
cone QCD sum rules. 



A Appendix 

In this appendix, the derivations of the explicit forms of the functions, /j, and Jj, appearing 
in Eq. (17) are presented and the method to calculate the continuum subtractions is 
demonstrated. For this purpose, consider a general term of the correlation function : 

n(p, q)^ [ d^xe''^^ /' du I Pa,e^("^+«"^)f-^(a,)/(«)^^T^?^^^ (Al) 

J Jo J {^y—X^)'' 

where Dai = daida2da3d{l — ai — a2 — a3) After switching to Euclidean space and carrying 
out the X integration, one obtains: 

U{p,q) = -i^ [du [ Vaidtt^-''ip{ai)f{u)e~^b^^'^'^'^ (A2) 

where — {q + kpf — q^k — p^kk + k{q + pY and k — a2 -\- was, k = I — k, and 
all the appearing momenta are Euclidean. In order to obtain Eq. (A2), the following 
representation of the Bessel function Ky[x) is used: 



Ky{mi,^x\) 1 foo dt \ rUb f x% 



1 /-"o di 

2 Jo V-^ 



/™2 \iy V In f'^+l 



exp 



(A3) 



where xe is the Euclidean position vector. Carrying out the Borel transformation with 
respect to {p + q)"^ using the identity: 

Se-"(-+^)^ = 5 - a) (A4) 

where B stands for the Borel transformation and is the Borel parameter and carrying 
out the t integral, one obtains: 



777,6 

where 



s{k) = <-'i'f^ + P'f^'k (A6) 



6 



where we have switched back to Minkowskian spacetime. Note that one can write the 
factors 1^ is the integrand as a differential operator acting on the exponential in the 
integrand and hence: 

n(M^) = -1^ (-^V [ duVa^e-^'^^^^^^ (AT) 
rub \ arrib J J k 

Note that 

e = / dse M^6(s — s(k)) (A8) 

which is nothing but the spectral representation of the exponential. The contributions of 
the higher states and the continuum are subtracted by replacing a cutoff, Sq, instead of the 
infinity as the upper limit of the spectral integral. Hence the effect of subtraction of the 
continuum and the higher states is to restrict the integration region in Eq. (A7) to the 
regions of k for which < s{k) < sq which implies that 6 < k <1 where 



ml + q^ -So + J {ml + - SqY - \ml{q^ -ml) 
5 = — ^- ^— (A9) 

where = mf^. Hence the contributions of the term given in Eq. (Al) to the sum rules 
after Borel transformation and the continuum is subtracted is given by: 

^/,.9x / d \^~'' f ^ ^ ip(ai)f(u) "(fc) , 

"<^^^ = -^i-ft;^j j iuVa,^-^^e-^ (AlO) 

where the integration region is determined by the conditions: 

8<k<l 

< w, ttj < 1 (All) 

and the Dirac's delta function in the definition of Pcuj fixes Z^j CKj = 1. 

In order to obtain the contributions to the sum rules from terms containing additional 
powers of one can use the trick of replacing by differentiation with respect to q^. For 
terms containing a factor of px in the denominator, one can use the following trick: in order 
not to have any singularity at = 0, the integral of these wave functions in the absence 
of the exponential should cancel. Hence, for these terms only, one can write: 



= ipx r dke^'^P'' (A12) 
Jo 



and the rest of the calculation is similar to the presented one. Note that the subtracted 1 
does not contribute. 

When the contributions of the continuum is subtracted and the Borel transformation 
is carried out for the other terms, the following functions which appear in Eq. (17) are 
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encountered: 



d 



druh 



rl — 6 pi — S — Qi 

/ dai / das 
Jo Jo 



rl rl—S rl—ai r 

du+ dai / / 

Jo Jo Jl-5-ai J- 



09(f\i. 1 — Hi — 03) «(l-«l-"3+"«3) 

■ -e 



1 — «! — as + ua^ 



d 



dm, 

-1-5 



1-5 
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1— 5— ai 



da- 



C du [ 
Jo Js 



l—al—as+uas 



dk 



pl—S pl—ai pi 

+ / dai / das L 

Jo Jl-S-ai J — 



du / dk 

"'a 



ip(ai,l — ai — as,as) »(i-«i-"3+""3) 



1 — «! — Qfs + uas 



d 



dm, 

-1-5 



1-5 







dai 



1— 5— ai 



da- 







du [ 

Jo Js 



1— al— aa+MQa 



dk 



+ 



rl — S rl — ai rl 

/ dai / das L 

Jo Jl-5-ai J — 



/•I— ai— as+uaa 

du / dk 

3 

_ 3(l-Qi-Q3 + tia3) 

(p{ai,l — ai — as,as)e 
du e 



dmh^ 
d ' 



u 



dm, 
d 



dk — : — e "m^ 



dmi, 



t du I 
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Figure 1: The dependence of the form factor f+{q ) on the Borel parameter M 
5, and 10 Gey 2 for so = 35 GeV^ and sq = 40 GeV^. 
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Figure 2: The same as Fig. (1) but for the form factor 
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Figure 3: The dependence of the form factor f+{(f') on at the continuum threshold 
So = 35 GeV^ and Sq = 40 GeV^ and at the Borel mass = 12 GeV^ and = 16 GeF^ 
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Figure 4: The same as Fig. (3) but for the form factor f-{q^). 



10 



References 

[1] Particle Data Group, C. Caso et. al Eur. Phys. J. C 15 (2000) 1 

[2] CLEO Collaboration: J. P. Alexander et. al, Phys Rev D61 (2000) 052001 

[3] P. Colengelo and A. Khodjamirian, "At Frontier of Particle Physics/ Handbook of 
QCD" ed. by M. Shifman (World Scientific, Singapore 2001), vol. 3, 1495 

[4] V. M. Braun, preprint hep-ph/9801222 (1998) 

[5] V. M. Braun, preprint hep-ph/9911206 (1999) 

[6] P. Colangelo and F. De Fazio, Phys. Lett B520 (2001) 78 

[7] C. S. Kim and Yo. Dong Yang, Phys. Rev. D65 (2002) 017501 

[8] V. M. Belyaev, V. M. Braun, A. Khodjamirian and R. Riickl, Phys. Rev. D51 (1995) 
6177 

[9] Tao Huang, Zho-Hong Li, Xiang-Yao Wu, Phys. Rev. D63 (2001) 1094001 
[10] F. De Fazio and M. R. Pennington JHEP 0007 (2000) 051 
[11] P. Ball, JHEP 9901 (1999) 010 

[12] T. M. Aliev,A. 6zpineci,M. Savci, Nucl. Phys. A 678 ; 034001(2001) , Phys. Rev D62 
:053012 (2000) 

[13] T. Feldmann, Int. J. Mod. Phys. A15, 159(2000) 

[14] T. M. Aliev and V. L. Eletsky Sov. J. Nucl. Phys. 38:936,1983 



11 



